Stabilization of the output power of intracavity frequency-doubled lasers 
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Intracavity frequency-doubled solid-state lasers exhibit intensity fluctuations of their light output, which are 
cause by nonlinear dynamical processes. Up to now, there are different solutions to this problem, but they reduce 
the output power, increase the size of the laser and/or make them more complicated to assemble. One focus of 
current research in nonlinear dynamics is derivation of control strategies from mathematical models and their 
experimental realization. We suggest a method to stabilize the output power by means of an electronic feedback 
of the emitted infrared light intensity to the pump power. First we show the theoretical predictions of a recently 
published stability analysis of a rate equation model with feedback. The presented experimental observation 
show systematic deviations from theory. This makes it necessary to refine the model to explain the deviations. 
The refinement has direct impact on the improvement of the feedback loop and, therefore, on the application of 
such a control scheme. 

PACS numbers: 05.45.Gg,42.55.Rz,42.55.Xi,42.60.Lh,42.60.Mi 
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Intracavity frequency-doubled solid-state lasers as light 
sources in the visible spectral range have a high technical po- 
tential because of their efficiency as well as their favorable 
ratio of compactness and output power QJ. Indeed, the intra- 
cavity frequency conversion design makes the nonlinear pro- 
cesses highly efficient, however, under multimode operation 
disadvantageous dynamical instabilities of the system, com- 
monly known as the green problem, are often observed in 
such kind of lasers. The instabilities originate from a cou- 
pling of the individual resonator modes due to the processes 
of sum frequency generation inside the nonlinear crystal and 
spatial hole burning within the laser-active media. Any kind 
of fluctuation of the output intensity can be a serious obstacle 
to many types of technical applications. Therefore, many peo- 
ple in industry and science are working on the stabilization of 
such kind of lasers. 

Since the first scientific discussion of the green problem by 
Baer |2] a number of proposals were made concerning the 
suppression of the fluctuations by means of different optical 
modifications of the resonator (e.g. ring resonator, mode se- 
lecting etalon)J3|,|4j,|5|]- An alternative approach was proposed 
by Roy et al. |6]. Leaving the optical setup of the system un- 
changed they exploited the existence of unstable steady states 
to perform a stabilization of the laser output intensity on a 
periodic orbit utilizing the so called occasional proportional 
feedback, known as OPF. As far as the application is con- 
cerned, it is more desirable to achieve a constant output in- 
tensity instead of a stable oscillating one. Nevertheless, it was 
demonstrated that it is possible to control the output intensity 
of the laser by a pure electronic feedback in principle. 

In this letter theoretical aspects and experimental results of 
an alternative control method are presented. From a rate equa- 
tion model 0| we derived a control scheme 0] that is able to 
stabilize the laser output by feeding back a control signal u(t). 
The signal u(t) can be derived from the sum intensities of the 
two orthogonal directions of the linear-polarized infrared light 



intensities which are emitted by the laser. More precisely, u(t) 
has the form 



U(t)=k x ( S ° x - Sx (t))+k y (S y -S y (t)), 



(1) 



where the light intensities s x and s y correspond to the sum of 
the individual modes polarized in the orthogonal directions x 
and y, respectively. The directions x and y correspond to the 
polarization eigenvectors of the laser cavity and depend on 
the birefringence and the length of the crystals inside the laser 
cavity, and are the time-independent sum intensities at 
the steady state and k x and k y are two control parameters. First 
experimental results on such kind of stabilization have been 
reported in 1 8 ] . 

The approach to form the feedback signal like Eq.Q re- 
sults from an analysis of rate equations, which consists of two 
ODEs for each individual lasing mode. Each set of ODEs 
for the individual laser modes incorporates a set parameter 
which are in general different from each other. The parame- 
ters are experimentally not accessible, therefore, for simplic- 
ity (w.l.o.g.) and to treat the system analytically we have cho- 
sen for all the laser modes the same set parameters. It can 
be shown, that the equations modelling the dynamics of the 
sum intensities s x and s y break off from the total system of 
equations describing the dynamics of each individual mode 
intensities in vicinity of the fixed point |7]. The stability of 
each individual laser mode intensity is reduced to the prob- 
lem of the stability of the sum intensities s x and s y . By using 
a control signal u(t) in accordance with Eq.Q and choosing 
proper values for the control parameters k x and k y one can the- 
oretically achieve and maintain the stability of the fixed point 
intensities as shown by Pyragas et al. |7J]. In the following we 
will give evidence that this control strategy also works exper- 
imentally in agreement with the theoretical predictions. 

The experimental setup of our intracavity frequency dou- 
bled Nd:YAG laser is shown schematically in Fig.^ A colli- 
mated beam of an infrared pump laser diode is focused on one 
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FIG. 1: Schematic drawing of the intracavity frequency doubled 
Nd:YAG laser pumped by a laser diode. 



end of a cylindrical Nd:YAG crystal. The optical resonator 
is formed by the planar HR coated surface of the Nd:YAG- 
crystal and a concave output mirror. A type II phase-matched 
KTP crystal is placed inside the cavity and serves as the opti- 
cal nonlinear element for generation of radiation at the second 
harmonics at A, =532nm. Due to the optical density of the 
Nd:YAG and the KTP crystal of about 1.82 and 1.8, respec- 
tively, the optical length of the resonator is about 43mm. 

The infrared and green spectral components of the total 
output intensity S was separated by a spectral selective beam 
splitter and the sum intensities s x and s y at 1064nm were sep- 
arated by polarizing beamsplitter. Special care was taken to 
operate the laser in the TEMoo mode of the laser resonator. In 
general, with increasing pump power the number of oscillat- 
ing longitudinal resonator modes increases. The dependence 
of the intensities 4 of the individual modes on the pump power 
Pio 1S complicated and is mainly governed by the process of 
spatial hole burning (J. Spatial hole burning originates from 
competition of modes for the same locally exited states inside 
the active media and results for instance in the fading of a laser 
mode with increasing pump power in favor of an other one. In 
order to decide which kind of nonlinear mode coupling takes 
place inside the KTP crystal, i.e. sum frequency generation 
or frequency doubling, the mode spectrum at A, =532nm can 
serve as a good indicator. The ratio between sum frequency 
generation and frequency doubling can be controlled by ro- 
tating the main axes of the birefringent KTP crystal |4]. For 
the following experiments a constellation with only sum fre- 
quency generation was chosen. 

Two photodetectors were used in order to detect the tem- 
poral changes of the sum intensities of the modes polarized 
in the two directions. The voltage signals of the detectors 
saved in a computer by means of an analog-to-digital con- 
verter (lOOMSa, 12bit). Two typical time traces are shown 
in the diagrams (a) and (b) of Fig. |2] The fluctuations dis- 
played in Fig. |3 a) are regular and were observed at a relative 
pump power of w = 1.89. The mode configuration was [7,4] 
(one mode polarized in x direction and four modes polarized 
in y direction). For the measurement in Fig. 0b) the relative 
pump power w was increased to a value of 2.07. All other 
system parameters remained unchanged. The corresponding 
mode configuration has changed to [1,5]. Irregular fluctua- 
tions can be seen. Note that the time traces of the sum inten- 
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FIG. 2: Measured time series of the infrared sum intensities s x , bold 
line, and s y , light line, polarized in x- and ^-direction, respectively, at 
two different relative pump powers w ((a) w = 1 .89 and (b) w = 2.07). 
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FIG. 3: Schematic diagram of the control feedback loop, to perform 
the stabilization by proportional feedback 



sity s x correspond to the dynamical behavior of a single mode 
in both diagrams and that the oscillations of s x and s y are out 
of phase indicating a mode competition process, known as the 
antiphase dynamics. For increasing pump power, the char- 
acteristic frequencies of the dynamics increase, following a 
square root behavior. The characteristic frequency is below 
200kHz for all the measurements presented here. 

The electric circuit used to implement the proportional 
feedback control is displayed schematically in Fig. [3] Instead 
of calculating the actual difference s* — s® (where * has to 
be x or y) as written in Eq. \l\ which is complicated, because 
in the experiment the values of the fixed point intensities 
are not known, we use a trick to estimate the difference by 
using high pass filtered signal. This approximation is valid 
because the fixed point intensity is equal to the time average 
of the intensity fluctuations. This procedure is discussed in 
more detail in |7J]. For this purpose we have used two high- 
pass filters with a cutoff frequency of about 1kHz. Subse- 
quently, the two signals are weighted with a variable gain am- 
plifier (— 1 < K xy < 1) and afterwards added up. The gain 
factors g x and g y of the amplifiers were adjusted by means 
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FIG. 4: Theoretical domains of stability (dark region) in the (k x ,k y ) 
plane obtained for proportional feedback control, (a) At a fixed mode 
configuration of [1,2] and a relative pump power of w = 1.38. (b) 
The mode configuration was [1,4] and the relative pump power was 
w = 2.07. 



of digital potentiometers and controlled remotely via a per- 
sonal computer. A third amplifier with variable gain K was 
used. Overall we end up with a voltage signal V c , given by 
V c = k x (V x - V/ c ) + k y (Vy - V* c ), where V* - V? c is approxi- 
mated by the high pass filtered signal from the detectors. Ad- 
ditionally, a limiter was used to restrict the distortions which 
were applied to the laser diode at high pump powers. This is 
only important during the time, where the control is switched 
on and the dynamics hasn't reached the vicinity of the fixed 
point. In order to feed the control voltage back into the laser, 
the current of the laser diode was modulated by a bypass con- 
nected in parallel to the laser diode. The -3dB bandwidth of 
the control circuit was about 300kHz and the phase shift was 
lower than 30 degrees at the highest characteristic frequency 
of the laser. 

In Fig.0]the theoretical domains of stability are displayed 
as shaded areas in the plane of the control parameters (k x ,k y ) 
according to theory 01 . The data were taken for two different 
exemplary situations with values of the relative pump power 
of w = 1.38 and w = 2.07 and mode configurations [7,2] and 
[1,4] in case of Fig.EJa) and Fig.01b), respectively. The pa- 
rameters have been chosen to be as close as possible to the 
performed experiments. Please note, that the mode config- 
uration is the only possibility in the model to have a handle 
on the mean sum intensity ration of the two polarization di- 
rections, because in the model all parameters for the different 
modes are chosen to be the same for simplicity. In the exper- 
iment it is more likely that the asymmetry of the intensities is 
caused by different sets of parameters for the different modes 
but they are experimentally not accessible. Most remarkable 
is the wedge-like shape of the two regions of stability of infi- 
nite size for larger k x and more negative k y values. Increasing 
the pump power, the angle between the two border lines of the 
region decreases. 

In Fig. 13 two measured domains of stability observed at 
different relative pump powers are displayed. We varied the 
two parameters (k x ,k y ) systematically by taking data at each 
combination of a hundred values in each direction and, there- 
fore, scanned the parameter plane. As a measure for the sta- 
bility, the standard deviation of the sum intensity s y was cho- 
sen, which directly relates to the amplitude of the fluctuations. 
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FIG. 5: Experimental domains of stability in the (k x ,k y ) parame- 
ter plane obtained for proportional feedback control (a) at a relative 
pump power w = 1.38 and the mode configuration [1,1] as well as (b) 
at a relative pump power w = 2.07 and the mode configuration [1,4]. 
The color codes the standard deviation of the laser intensity output 
(the darker the color, the lower the standard deviation, i.e., black is 
closed to zero deviation). The two straight bold lines in both of the 
diagrams are added to guide the eye. 



Darker regions correspond to lower values of the standard de- 
viation. Thus, the black regions mark the domain with nearly 
no fluctuations, that is, the domain in which the laser has a 
constant output power. For the two measurements shown in 
Fig- HI only the pump power was changed, all other system 
parameters were kept constant. For the measurement corre- 
sponding to Fig. |2a) the relative pump power w was 1.38. 
The laser oscillated in two orthogonal modes, [7,7] and the 
output intensity showed weak periodic fluctuations. For the 
second measurement corresponding to Fig. |5jb), the pump 
power was increased to a value of w = 2.07. At this pump 
power the fluctuations were highly irregular with a dynamic 
range of 100%. This means, the laser switched on and off. 
The mode configuration was [1,4]. As one can see in Fig. [5] 
the output power could be stabilized at different pump pow- 
ers and different mode configurations. The accordance of the 
wedge-like shape of the experimental and the theoretical do- 
mains of stability is interesting. Also the tendency of the angle 
of the wedge (indicated in Fig.[5]by two black lines) to de- 
crease with increasing pump power is nicely reproduced. But 
in contrast to the theoretical predictions, the domains of sta- 
bility are finite and shrink in general due to increasing pump 
power to zero (in the example presented here, the stabilization 
fails entirely for w > 2.2). An extension of the fundamental 
model didn't give the right explanation fioll . 

To understand this behavior we have to consider an impor- 
tant issue, which we have neglected up to now. The region 
of stability was calculated using linear stability analysis by 
only looking at the real part of the eigenvalue without taking 
into consideration the imaginary part of the eigenvalue. Look- 
ing at the imaginary part of the eigenvalue we notice that the 
imaginary part increases if k y becomes more and more nega- 
tive. From the experimental point of view it is obvious that the 
feedback signal is not modelled correctly. In any experimen- 
tal situation the bandwidth of the feedback loop is restricted 
which was not considered in the model — in the model the 
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FIG. 6: The domains of stability in the (k x ,ky) parameter plane ob- 
tained from the analysis of the model for proportional feedback con- 
trol with a restricted bandwidth at (a) rel. pump power w = 1.38 at 
a mode configuration [1,2] and (b) rel. pump power w = 2.07 at a 
mode configuration [1,4], The cutoff frequency of the low-pass filter 
was chosen similar to the experimental one to be 300kHz. The lines 
are marking the position of the boundaries of the domains of stability 
for an infinite bandwidth. 



bandwidth is unlimited. This fact influences not only the re- 
gion of stability in the (k x ,k y ) parameter plane, but also the 
behavior of the stability region under variation of the pump 
power. As we already mentioned above, the characteristic fre- 
quency of the fluctuations of the laser output intensity grows 
by increasing the pump power and, therefore, the region of 
stability also shrinks for a certain bandwidth of the feedback 
loop if the pump power is increased. 

For the stability of the system it is most important to take 
an upper cutoff frequency into account. The model has to be 
modified by an additional degree of freedom introduced by 
the low pass characteristics of the feedback loop m/. This can 
be modelled by adding an additional differential equation to 
the rate equations introduced in |7]. The additional equation 
reads as follows: 



dm 



*('))• 



M4- 



-Sy(t))-Ul), 



(2) 



where Cfl c denotes the cutoff frequency of the low pass filter. 
Using this form of the feedback signal and performing the sta- 
bility analysis we obtain the results presented in Fig. |6] For 
the simulation a cutoff frequency of 300kHz was used, which 
corresponds to the experimental situation. It can be seen that 
now all essential features of the experimentally observed re- 
gion of stability can now be reproduced. This finding shows 
that an increased bandwidth gives the possibility to stabilize 
the laser even further to higher pump powers up to values im- 
portant for technical applications. 



In conclusion we presented experimental evidence that an 
intracavity frequency-doubled solid state laser can be stabi- 
lized by a proportional feedback of the two sum intensities of 
the orthogonal polarized modes to the pump power. By refin- 
ing an earlier model it was possible to understand the shape 
and dependence on system parameters of the region of stabil- 
ity in the control parameter plane. The presented result aim 
in the direction what to do for an application of the method 
in a laser with more realistic values of the pump power and, 
therefore, of the intensity of the light output. To extend the 
applicability of the method to a hundred times higher pump 
powers the bandwidth has to be increased by the factor ten. 
This can be easily achieved. We are convinced that these find- 
ings are of great importance for technical application of such 
a kind of laser. 
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